High-dielectric constant and wide band gap inverse silver oxide 
phases of the ordered ternary alloys of Si02, Ge02 and Sn02 
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High-dielectric constant and wide band gap oxides have important technological applications. 
The crystalline oxide polymorphs having lattice constant compatibility to silicon are particularly 
desirable. One recently reported candidate is the inverse silver oxide phase of Si02. First-principles 
study of this system together with its isovalent equivalents Ge02, Sn02 as well as their ternary 
alloys are performed. Within the framework of density functional theory both generalized gradient 
approximation and local density approximation (LDA) are employed to obtain their structural 
properties, elastic constants and the electronic band structures. To check the stability of these 
materials, phonon dispersion curves are computed which indicate that Ge02 and Sn02 have negative 
phonon branches whereas their ternary alloys Sio.5Geo.5O2, Sio.5Sno.5O2, and Geo.5Sno.5O2 are all 
stable within LDA possessing dielectric constants ranging between 10 to 20. Furthermore, the 
lattice constant of Sio. 5660.502 is virtually identical to the Si(lOO) surface. The GW band gaps of 
the stable materials are computed which restore the wide band gap values in addition to their high 
dielectric constants. 

PACS numbers: 61.50.Ah, 62.20.Dc, 63.20.Dj, 71. 15. Mb, 71.20.-b, 71.20.Ps, 77.22.Ch 



o 

> 
m 
vn 
m 
i> 
o 
\o 
o 



I 

o 
o 



High-dielectric constant and wide band gap oxides are 
of general interest for the next-generation gate oxides for 
silicon-based electronics^ and also as host matrices for 
nonvolatile flash memory applications.^ Amorphous ox- 
ides have been generally preferred as they are good glass- 
formers which tend to minimize the number of dangling 
bonds at the interface. In this respect, poly-crystalline 
oxides are undesirable as the grain boundaries cause 
higher leakage currents and possible diffusion paths for 
dopantsii On the other hand, crystalline oxide grown epi- 
taxially on silicon^ can be favorable as it will result in 
high interface quality provided that it is lattice-matched 
to Si. 

Very recently, Ouyang and Ching* have reported a 
high-density cubic polymorph of Si02 in the inverse 
Ag20 structure, named by them as the i-phase, possess- 
ing both high dielectric constant, as in stishovite phase, 
and the lattice constant compatibility to Si(lOO) face 
which make it very attractive for electronic applications. 
In this computational work, we continue this search for 
the crystalline high-dielectric constant oxides with the i- 
phases of Ge02 and Sn02 as well as their ordered ternary 
alloys with Si02. This pursuit is in line with the Inter- 
national Technology Roadmap for Semiconductors where 
computational synthesis of novel high-dielectric materials 
is emphasizedj^ We employ the well-established ah initio 
framework based on the density functional theory within 
the generalized gradient approximation and local density 
approximation using pseudopotentials and a plane wave 
basis 

The unit cell for the ordered ternary alloy X0.5Y0.5O2 
in the inverse Ag20 structure is shown in Fig. [T] Struc- 
tural and electronic properties of the i-phase structures 
under consideration have been calculated within the 
density functional theory;^ using the plane wave basis 
pseudopotential method as implemented in the ABINIT 



code.'' The results are obtained under the the general- 
ized gradient approximation (GGA) and local density ap- 
proximation (LDA) where for the exchange-correlation 
interactions we use the Teter-Pade parameterization)^ 
which reproduces Perdew-Zunger^ (which reproduces the 
quantum Monte Carlo electron gas data of Ceperley and 
Aider™). We tested the LDA results under two differ- 
ent norm-conserving TrouUier and Martinsii type pseu- 
dopotentials, which were generated by A. Khein and D.C. 
AUan (KA) and Fritz Haber Institute (FHI); for either 
set, the d electrons were not included in the valence con- 
figuration. Our calculated values for these two types of 
pseudopotentials were very similar. 




FIG. 1: Ball and stick model of the i-phase ordered ternary 
alloy X0.5Y0.5O2. 

In the course of both GGA and LDA computations, 
the plane wave energy cutoff and k-point sampling were 
chosen to assure a 0.001 eV energy convergence for all 
i-phase crystals. In the case of Si02 this demands a 
65 Ha plane wave energy cutoff and 10x10x10 fc-point 
sampling. Phonon dispersions and phonon density of 
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states were computed by the PHON prograrc^^ using 
a 2 X 2 X 2 supercell of 48 atoms to construct the dy- 
namical matrix. The required forces were extracted from 
ABINIT. The corrected band gap values are computed by 
obtaining self-energy corrections to the DFT Kohn-Sham 
eigenvalues in the GW approximation.^'^ All parameters 
used during the GW calculation were chosen to assure a 
0.001 eV energy convergence. 



TABLE I: First-principles LDA and GGA structural data for 
i-phase crystals. 

Crystal a (A) Density (gr/cm^) x-O (A) y-O (A) 



SiOa 


LDA 3.734 


3.830 


1.617 






GGA 3.801 


3.633 


1.646 




GeOa 


LDA 3.916 


5.781 


1.696 






GGA 4.053 


5.215 


1.755 




Sn02 


LDA 4.180 


6.864 


1.808 






GGA 4.452 


5.671 


1.928 




Geo.5Sio.5O2 


LDA 3.836 


4.843 


1.697 


1.625 




GGA 3.923 


4.528 


1.762 


1.635 


Geo.5Sno.5O2 


LDA 4.042 


6.416 


1.688 


1.813 




GGA 4.250 


5.522 


1.748 


1.932 


Sno. 5810.502 


LDA 3.970 


5.590 


1.818 


1.620 




GGA 4.114 


5.015 


1.935 


1.628 



Using XO2 and X0.5Y0.5O2 as the generic notation, 
the 0-X-O and 0-Y-O bond angles are 109.47° and the 
X-O-X and X-O-Y bond angles are 180° according to 
the crystal construction of this cubic i-phase (cf. Fig. [1]). 
Other structural information such as the lattice constants 
and bond lengths of all i-phase crystals are listed in Ta- 
ble m The Si(lOO) surface lattice constant is about 



TABLE II: Elastic constants and bulk modulus for each crys- 
tal. 



Crystal 




Cii(GPa) 


Ci2(GPa) 


C44(GPa 


B(GPa) 


Si02 


LDA 


383.6 


260.0 


243.0 


301 




GGA 


354.3 


232.1 


227.9 


273 


Ge02 


LDA 


297.0 


231.2 


175.6 


253 


Sn02 


LDA 


208.9 


185.5 


113.9 


193 


Geo.5Sio.5O2 


LDA 


349.4 


253.2 


200.0 


285 




GGA 


292.8 


203.9 


161.8 


234 


Geo.5Sno.5O2 


LDA 


255.4 


210.8 


106.3 


226 


Sno. 5810.502 


LDA 


277.5 


217.4 


103.9 


237 




GGA 


238.3 


183.0 


202.8 


201 



3.83 A, therefore according to LDA results Sio.5Geo.5O2 
is of particular interest as it can be epitaxially grown on 
Si without any strain. According to our well-converged 
calculations Sio.5Geo.5O2 has a lower total energy com- 
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(a) SiO, 




DOS 

FIG. 3: LDA electronic band structure and DOS (States/eV 
(d) Sno.5Sio.5O2. 

pared to both Si02 and Ge02, the latter itself is unstable 
as will be shown later; this can be taken as some indica- 
tion of immunity to the phase separation of this ternary 
alloy into its binary compounds. 

The LDA and GGA results of the three independent 
elastic constants and bulk modulus for all crystals are 
tabulated in Table HIl An important concern is the sta- 
bility of these cubic phases. The requirement of mechan- 
ical stability on the elastic constants in a cubic crystal 
leads to the following constraints: Cn > C12, Cn > 0, 




DOS 

ill) of i-phase (a) Si02, (b) Geo.5Sio.5O2, (c) Geo.5Sno.5O2, and 



C44 > 0, and Cii + 2Ci2 > 0. The clastic constants cal- 
culated by both LDA and GGA shown in Table |TT] sat- 
isfy these stability conditions. Furthermore, we compute 
the LDA and GGA phonon dispersion curves of these 
structures using the PHON program. First, to verify 
the validity of the results of the PHON program we com- 
pute the phonon dispersions of the Si02 and Ge02 by us- 
ing both PHON and ANADDB extension of the ABINIT 
code. There exists a good agreement between two calcu- 
lations. Next, we calculate the phonon dispersions of the 



4 



all i-phase crystals via PHON program with forces ob- 
tained from LDA and GGA. It is observed that Si02 is 
at least locally stable whereas Ge02 and Sn02 contains 
negative phonon branches which signal an instability of 
these phases. As for their alloy, Geo.5Sno.5O2, accord- 
ing to LDA this material is stable whereas within GGA 
it comes out as unstable. For the stable structures the 
LDA phonon dispersions and the associated phonon den- 
sity of states (DOS) are shown in Fig. [2] 



TABLE III: LDA and GGA dielectric permittivity tensor for 
the stable crystals. 



Crystal 




^xx yy zz 


00 00 00 

^xx yy zz 


SiOa 


LDA 


9.857 


3.285 




GGA 


9.970 


3.303 


Geo.5Sio.5O2 


LDA 


11.730 


3.416 




GGA 


14.383 


3.585 


Geo.5Sno.5O2 


LDA 


19.415 


3.527 


Sno.5Sio.5O2 


LDA 


12.883 


3.360 




GGA 


18.096 


3.711 



For the stable systems, the static and high-frequency 
dielectric constants are listed in Table IIIII The static 
dielectric constants falling in the range between 10 to 
20 suggest that these are moderately high dielectric con- 
stant crystals. It can be observed that GGA yields sys- 
tematically higher values for the dielectric constants of 
these structures. Employing KA pseudopotentials, the 
LDA band structure for the crystals are displayed along 
the high-symmetry lines in Fig. [3] including the elec- 
tronic DOS. The widths of the valence bands get progres- 
sively narrowed from Fig. [31[a) to (d), i.e., from Si02 to 
Sno.5Sio.5O2. For all of the i-phase crystals under consid- 
eration including the unstable ones the conduction band 
minima occur at the F point whereas the valence band 
maxima are located at i? point making them indirect 
band gap semiconductors. As tabulated in Table HVl the 
direct band gap values are only marginally above the in- 
direct band gap values. Again GGA systematically yields 
narrower band gaps compared to LDA. 

A renown artifact of LDA is that for semiconductors 
and insulators band gaps are underestimated.'' In this 
work, the corrected band gap values are also provided by 
GW approximation. As there are different GW imple- 
mentations we briefly highlight the particular method- 
ology followed in the ABINIT code. First, a converged 
ground state calculation (at fixed lattice parameters and 
atomic positions) is done to get self-consistent density 
and potential, and Kohn-Sham eigenvalues and eigen- 
functions at the relevant band extrema /c-points as well as 
on a regular grid of fc-points. Next, on the basis of these 
available Kohn-Sham data, the independent-particle sus- 
ceptibility matrix xo is computed on a regular grid of 
g-points, for at least two frequencies (usually, zero fre- 
quency and a large pure imaginary frequency - on the 
order of the plasmon frequency, a dozen of eV). Finally, 



the Random Phase Approximation susceptibility matrix, 
X, the dielectric matrix e and its inverse are com- 
puted. On this basis, the self-energy, S matrix element 
at the given fc-point is computed to derive the GW eigen- 
values for the target states at this fc-point. Note that this 
GW correction is achieved as a one-shot calculation (i.e., 
no overall self-consistency) hence, our results technically 
corresponds to Go Wo which has been the standard ap- 
proach as originally proposal by Hedin.^^ The GW cor- 
rection as can be observed from Table IIVI restores the 
wide band gap values; this feature is essential for these 
materials to provide sufficient confinement to carriers of 
the narrow band gap semiconductors such as silicon. 



TABLE IV: Indirect {Eg) and direct {Egx) band gaps for 
each i-phase crystal within LDA, GGA, and for the stable 
structures the GW approximation (GWA). 



Crystal 




Eg (eV) 


Eg.y (eV) 


Si02 


LDA 


5.269 


5.870 




GGA 


4.584 


5.155 




GWA 


7.283 


7.964 


Ge02 


LDA 


2.402 


2.511 


Sn02 


LDA 


2.285 


2.670 


Geo.5Sio.5O2 


LDA 


3.666 


4.179 




GGA 


2.558 


3.005 




GWA 


5.943 


6.513 


Geo.5Sno.5O2 


LDA 


2.487 


2.900 




GGA 


0.767 


0.865 




GWA 


4.533 


4.972 


Sno.5Sio.5O2 


LDA 


3.292 


3.900 




GGA 


1.763 


2.304 




GWA 


5.484 


6.153 



We have also considered the i-phase of Pb02 which 
turned out to be unstable and hence its ah initio data are 
not included. In this work, we do not consider the ther- 
modynamic stability of these i-phase oxides. However, 
for technological applications rather than bulk systems 
the epitaxial growth conditions become more critical.^"' A 
promising direction can be the finite temperature inves- 
tigatioi>i^ of these i-phase isovalent structures on Si(lOO) 
surfaces using large number of monolayers. 

This first-principles study suggests that the i-phases of 
Ge02 and Sn02 are unstable whereas Si02, Sio.5Geo.5O2, 
Sio.5Sno.5O2 are particularly promising due to their high 
dielectric constants as well as wide band gaps as re- 
stored by the GW correction. Moreover, they are 
lattice-matched to Si(lOO) face, especially for the case of 
Sio.5Geo.5O2. We beheve that these findings can further 
boost the research on the crystalline oxides. 
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